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OVER CONVECTION- AND TRANSEERATION-COOLED SURFACES 


By John N. B. Livlngood and Patrick L. Donoughe 


SUMMARY 

A summary of exact solutions of the laminar-houndary-layer equations 
for vredge-type flow, useful in estimating heat transfer to such arbitrar- 
ily shaped bodies as turbine blades, is presented. The solutions are de- 
termined for small Mach numbers and a Prandtl number at the wall of 0.7j 
ranges of mainstream pressure gradients and rates of coolant flow through 
a porous wall are considered for the following cases: (l) small tempera- 

ture changes in the boundary layer along a constant- and along a vaidable- 
temperature wall, and (2) large tenperature changes in the boundary layer 
along a constant -temperature wall. 

Dimensionless forms of heat-transfer and friction parameters and 
boundary -layer thicknesses are tabulated. Q3ie results indicate that cool- 
ant emission and increased stream-to-wall tenperature ratios diminished 
the friction and heat transfer for a constant wall temperature. For a 
variable wall temperat\jre with small temperature differences in the bound- 
ary layer, the friction was unaffected, but the heat transfer was greatly 
increased for increased wall-temperature gradient. Heat-transfer resialts 
in the literature reveeil that transpiration cooling is much more effective 
for Prandtl numbers of the order of 5.0 than for 0.7. 


INTRODUCTION 

The cooling of structural parts in propulsion systems such as gas- 
tiorbine blades, combustion-chamber walls, and rocket nozzles has become 
increasingly important with the growing demand for hi^er-powered and more 
efficient systems. Such cooling may be accomplished by either convection 
or transpiration methods. Convection cooling is obtained by passing cool- 
ing air along the coolant side of the wall; transpiration cooling is ac- 
complished by forcing the coolant throu^ a permeable wall. Methods pre- 
viously used to predict heat transfer had to be extended to include 
transpiration cooling for the different types of application. The NAGA 
Lewis laboratory has Initiated a program to extend such methods in the 
laminar -flow region. 
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Exact solutions of laminar -boundary-layer equations for wedge-type 
flow have been found usefxil in predicting laminar heat transfer to bodies 
with impermeable walls (ref. l) . (Wedge-type flow is flow for which the 
mainstream velocity is proportional to a povrer of the distance from the 
stagnation point.) However, only a few solutions for heat transfer ex- 
isted for transpiration-cooled wedges (ref. 2); these were determined for 
conditions of constant wall temperature, constant flvdd properties, and 
a Prandtl number of vinity. In order to apply this method to the cooling 
of sucb structural parts of propTilsion systems as turbine blades, wedge 
solutions were required for conditions applying to sucb components. These 
conditions include the slmultaneoTis effects of flow throu^ a porous wall, 
large pressirre gradients in the mainstream, small Mach niimbers, small tem- 
perature differences throu^ the laminar boundary layer (constant property 
values or a stream-to-wall temperatujre ratio near unity) , and large tem- 
perature differences through the lami nar boundary layer (variable proper- 
ty values or stream-to-wall tenperature ratios appreciably different from 
unity) . 

Exact solutions of the laminar-boundary-layer equations for a range 
of flows throu^ a porous wall, for a range of pressure gradients, and 
for several stream-to-wall absolute temperature ratios are presented in 
references 3 to 5 for a Prandtl niamber at the wall of 0.7. In these solu- 
tions the wall temperature was assumed constant, and the range of stream- 
to-wall tenperature ratios was chosen to include the cases both of con- 
stant property values and of property values that varied with powers of 
the absolute tenperature. Wedge solutions are presented in reference 6 
for constant property values and a variable wall temperature. Again the 
Prandtl number was chosen as 0.7. Solutions for isothermal impermeable 
flat plates are contained in inference 7 for various Prandtl numbers; 
the results are extended to wedge-type flow in reference 8. Solutions 
for inp>ermeable wedges that illustrate the effect of Prandtl number, pres- 
sure gradient, and surface temperature variation are presented and sum- 
marized in reference 9. Reference 10 presents analytical heat-transfer 
results for an isothermal permeable flat plate with fluids of different 
Prandtl numbers. 

Wedge solutions can be utilized to obtain approximate Information 
on heat transfer to bodies of arbitrary shape. A first approximation can 
be obtained by stipulating that the heat-transfer coefficient at any loca- 
tion on a body of arbitrary cross section is identical with the heat- 
transfer coefficient on a wedge for which, at the same distance from the 
stagnation point, the stream velocity and the stream- velocity gradient 
ara the same as those on the arbitrarily shaped body. (Such use of wedge 
solutions is nade for an inpermeable cylinder in ref. 1. The method is 
described for permeable bodies in appendix C of ref. 11.) This approxi- 
mation satisfies local straam conditions only, and does not properly ac- 
count for the development of the boundary layer ahead of the point under 
consideration. Taki.ng into account the development of the boundary layer 
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upstream of a point under consideration leads to a second and letter ap- 
proximation for heat transfer to a "body of arhitrary cross section. In 
this method it is assumed that the rate of increase of the houndaiy-layer 
thickness is the same at a location on the arhitrarily shaped "body as at 
the position on a wed^ that has the same boundary-layer thickness, stream 
velocity, and stream- velocity gradient. This method is proposed for im- 
permeable bodies in reference 8, is extended to transpiration cooling with 
small tempei^ture differences in reference 12, and to transpiration cool- 
ing with large ten 5 )erature differences in reference 11. Eeference 11 al- 
so contains a complete summary of. other types of laminar-boundary -layer 
solutions. The calculation of a transpiration-cooled turbine blade using 
the method of reference 11 is presented in reference 13. A review of some 
of the available laminar-boimdary-l^yer solutions is given in table I in 
chronological order. 

The two methods that employ wedge solutions for predicting heat trans- 
fer to bodies of arbitrary cross section are compared with available ex- 
perimental data in reference 11. For an impermeable circular cylinder, 
heat-transfer coefficients determined directly from wedge solutions agreed 
within 15 i>ercent ’t/ith the results of experiments. Use of the second ap- 
proximation inproved the agreement to 8 percent. Direct Tose of wedge 
solutions for predicting heat transfer to a symmetrical airfoil with an 
8-foot chord resulted in agreement within the experimental scatter of the 
data (ref. 14). 

From the foregoing discussion, it is apparent that laminar-boundary- 
layer solutions for wedge-type flow are useful in ceilculating heat trans- 
fer to bodies of arbitrary shape, such as cylinders, airfoils, turbine 
blades, and so forth. It is also apparent that the literature contains 
many investigations on wedge-type flow. Including results for fluids with 
different Prendtl numbers florfng through porous flat plates. It is the 
purpose of this report to collect and summarize ttiese wedge-type flow 
analyses for use in predicting heat transfer to bodies of arbitrary shape, 
and to compare heat-transfer results for fluids with different Prandtl 
numbers flowing throu^ porous flat plates. 


SYMBOLS 


The following symbols are used in this report: 


B,C constants of proportionality 


C.p „ ■ local skin-friction coefficient baaed on properties at wall 

X ,w ^ 

temperature, p 

PwU„ 
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00 


C 

p 


Eu 


l,ocal skLn-frictioii coefficieut "based on properties at stream 

temperature, 2 

P u 
00 00 

~~ 2 ~ 

specific heat at constant pressure 

-X ^ 

Euler nmher, —} U for wedge-type flow 

P 

00 00 


f . dimensionless stream function, 

f ‘,f",f"’ first, second, and third derivatives of f with respect to 

H local heat-transfer coefficient 

k thermal conductivity 

Eu loc^ Husselt number, Bx/k^ 


n 


Pr 

P 

Re 


X ^ 

waH-teraperatxire-greidient parameter, ' i" ' ^ 

5 !^ - 


Prandtl number, 
static pressure 

U 

local Reynolds number, 


T ten^erabure 

fluid velocity at edge of boundary layer 

u fluid velocity in boundary layer parallel to wall 

V fluid velocity in bomdary layer normal to wall 
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X 

y 

a 

5 

5* 


5 


c 


5 


i 


5 


t 


e 


n 

9 

e*,e" 

n 

V 

p 

T 


distance along surface 
distance normal to surface 

exponent of temperature for specific tieat, c^ oc aP’ 


■boundary-layer thickness 


displacement-layer thickness. 



convection- layer thickness. 


pu 


- V 


dy 


momentum- la 3 rer thickness. 



a 


thermal -layer thickness. 




T - T 

OO 

5^ ay 


exponent of temperature for thermal conductivity, k oc 


dimensionless boundary-layer coordinate; 




temperature -difference ratio. 




first and. second, derivatives of 9 with respect to ^ 
absolute viscosity of fluid, 
kinematic viscosity, p/p 
density of fluid 


shear stress 
stream function 
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05 esqjonent of ten^jeratirre for viscosity, p. (x 

Subscripts: 
w ■wall 

00 mainstream, outside 'boundary layer 


AMLYSES 

Solutions of the laminar-'boundary-layer equations for small -ten^jera- 
ture changes in the houndary layer along constant- and. variable - 
temperature walls and. for large temperature changes in the boundary layer 
along a constant-temperature wall are discussed.. Details of the analyses 
and C8G_culation methods asre not contained herein, but are available in 
refeiences 3 to 6. 


Boundary -Layer Equations 

The equations of the laminar "boundary layer for steady-state flow 
of a -viscous flioid -with heat transfer are as follows; 

Momentum: 

( 1 ) 

Continuity: 

^ (pu) + ^ (pv) = 0 (2) 

Energy: If the temperat'ure differences between the -wall and the 

mainstream are assumed large comrpared -with temperature changes caused by 
compression and frictional heating, the energy equation may be written 



Equations (l) to (3) include transpiration cooling for the case when 
the same medium as that in the outside flow is used as coolant and the 
following boundary conditions are defined: 


pu 




5y 


¥ 




Su\ 
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and 


u 


= 0, V = and T = for y = 0 


u -»• U and T -+• T for y -> «> 

00 00 ^ 


( 4 ) 


QJbe property values p, c ^ and p appearing in the equations are as- 
sumed to he either constant or functions of temperature only. 


Change of Variables 


A transformation of the partial differential equations is accom- 
plished hy changing the independent variables x and y into new in- 
dependent variables x and t] where 



( 5 ) 


In addition^ the dependent variables u, v^ and T are replaced by the 
new dependent variables f and 9 (refs. 15 and 7), where 


Pw^ 




f = 




T - T. 


e = 


w 

T - T 
•^00 -^w 


( 6 ) 




The continuity equation is satisfied by the stream function ijr, and hence 


pu = 






( 7 ) 


Application of these changes in variables to equations (l) to (3) 
results in a set of equations that, in principle, can be solved for two- 
dimensional flow about any arbitrary shape. For engineering piurposes, 
however, obtaining such solutions is not feasible. As a consequence, ex- 
act solutions for wedge-type flow will be presented herein that can be 
used by the two methods e^^lained in the XNTROnJCTION to approximate heat 
transfer to bodies of arbltreiry cross section. 
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Ordinary Differential Equations for Small 
Stream-to-W^ll Temperature Differences 
(Constant and. Variable WeOLl Temperatures) 

As already stated., for constant property values including density 
(small temperature differences), "both constant and. variable wall tempera- 
tures are considered, herein. For variable Trail temperatTure, it is as- 
sumed. that the difference between the wall and. stream temperatxires is 
proportional to a power of the distance from the leading edge^ that is, 

T^ - T„ = (8) 


where n and. T„ are constants (refs. 16 to 18).’ Differentiation of 

equation (8) yields the following expression for the wall-tenperature- 
gradient parameter: 


n 



dx' 


( 9 ) 


For wedge -type flair the stream velocity varies as a power function 
of the distance from the stagnation point measured, along the svirface: 

= Cx® (10) 


It is customary to refer to the exponent m in this equation as Euler 
number. The Euler number can be expressed, by the Bernoulli, equation 


dp 

„ ^ X ^oo 

= U dx 
p IP W 
00 00 

~X 

Assuming that f and 9 are functions of t] only and iising equa- 
tions (5) to (U) transform equations (l) to (3) into the total differen- 
tial equations; 



f t f I 


Euf 


Eu + 1 
2 


ffi. 


- Eu 


Eu + 1 


Prf0' 


nPrf>(l - 0) 


( 12 ) 


0 " 


2 


(13) 
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with the houndflry conditions 

f = f . f ’ = 0, and. 9 = 0 for ri = 0 

and 

f ' 1 and 0 -*• 1 for tj -+■ ,oo 

Equation (13) is identical to equation (l2) of reference 6 when the suh- 

T - T 
00 

stltution y = 1 - 0 is made, where Y = ^ ^ . 

W "^00 

The velocity at ttie surface (y = O) is 
w 

I 

Eu + 1 
^w 2 

The transformation, therefore, prescrihes a certain variation of the cool- 
ing velocity v^ along the surface, since the function f^ has to he 

constant (independent of x) . Because the stream velocity is described 
hy equation (lO), the coolant velocity is also proportional to some 

2 

power of X, namely oc x . Reference 2 shows that such a varia- 
tion of the coolant velocity leads to a constant wall temperature when 
heat transfer hy conduction and radiation may he neglected. Only conduc- 
tion along the wall, or radiation, or hoth may lead to the case of a wall- 
temperature variation if f.^ is constant. 


Ordinary Differential Equations for Large Stream- 
to-Wall Teraperattnre Differences 
(Constant Wall Ten^jerature) 

In order to account for changes in the fluid properties, the follow- 
ing property variations with absolute tenqperat\ire are employed: 

p oc 0^, k oc T^, oc of, p oc (l6) 

Use of equations (16) and equations (5), (6), (7), (lO), and (ll) trans- 
forms equations (l) and (3) into the total differential equations: 
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Eu + 1 

0" = - —2— 


f = Eu 


(I;) 

^-<0 r V - 

f^\ 2 Eu + 1 /^\ 

V^y 2 


/^„\ /m v-tO-1 

.1 "■-“(€) © 


Eu 


■ ‘)©"‘ "■•■ - ft - W 


f'O' 




'T, 


-2 


CD 


“ - ll (JL\ f'0'2 




The houndary conditions are 


© 


= o 


f = f'=0f and 0=0 for t] = 0 


and 




f ’ ~ and 0 -► 1 for t) » y 

CO 


The velocity at the surface is now given hy 


^w “ - 2 


Eu + 1 |^■^Poo 




(18) 


(19) 


Boundary-Layer Thicknesses, Heat Transfer, and Friction 
From equations (5) to (7) it follows that 

pu 


Pw^oc 


= f ’ 


The speciflc-wei^t-flow lutlo is tlierefore expressihle as 

T 

pu f I «> 
p U T 

^ rv-i TJ 


( 20 ) 
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and the velocity ratio is expressible as 





( 21 ) 


From equations (5) to (v), (20), and (21)^ and tbe definitions of 
the boundary-layer tbicknesses given in the SYMBOLS, dimensionless forms 
of the displacement, momentum, convection, and thermal boundary-layer 
thicknesses may be expressed as foUowa: 



( 22 ) 


(23) 


(24) 


(25) 


Equating the convective heat transfer H(T^ - T^) to the conductive 
heat transfer k^(3T/Sy)^ and using equations (5) and (6) yield the 
dimensionless form of the heat-transfer parameter 


Nu 

/Vse 


0 * 

w 


With the wall shear stress given by 


T = 

w 



(26) 


( 27 ) 


and friction coefficients defined as 
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or 




w 


P 

W 00 


T 

Poo^oo 

2 


the dimensionless forms of the friction parameter are 



or 



(28) 


Assumptions for Numerical Calculations 

The numerical solutions veve obtained, for a Erandtl number at the 
wall of 0.7 (appropidate for air), pressure variations represented by 
values of Eu of 0, l/2, and 1, flow rates throu^ the porous wall rep- 
resented by values of f^ of 0, -l/2, and -1, and wall-temperature varia- 
tions represented by values of n from the value corresponding to a zero 
wall-temperature gradient to unity. 

For large temperature differences, the exponents in equation (16) are 
required. These exponents, determined by logarithmic plots of the proper- 
ties of air against temperatTure (600° to 2400° F) izslng the property val- 
ues given in reference 19, are as follows: 


Exponent in viscosity -temperature relation, co 0.7 

Exponent in thermal-conductivity - ten 5 >eratiire relation, e .... 0.85 
Exponent in specific-heat - ten^erature relation, a 0.19 


EESULTS AND DISCUSSION 

The heat-transfer and friction parauKters and bo\mdary-layer thick- 
nesses are presented in table H. In addition, values of n that result 
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in = 0 (for small temperatTore differences) and values of Eu that 
result In f^" = 0 (for constant wall temperature) are also given. Val- 
ues of f , 0, and their derivatives may be obtained from references 4 to 
6 except for cases \diere f^ = 0, n = 0, and Eu < 0. For these cases, 

only f ' and 6 are listed. Representative profiles are shown In fig- 
ures 1 to 4 herein. Unless otherwise noted, the source of the distribu- 
tion Is given In table II. 

(O 

Boundary-Layer Profiles 

Flgixre 1 presents velocity distributions In the boundary layer. 

Parts (a), (c), and (e) are for the flat plate alined with the mainstream 
(Eu = 0) j parts (b) , (d) , and (f ) are for the plate perpendicular to the 
stream or stagnation flow (Eu = l) . Figures l(a) and (b) illustrate the 
effects of coolant emission or blowing throu^ the porous wall for con- 
stant property values. Distributions for additional blowing rates are 
obtained from references 20 and 21. 

- Because eguatlons (12) and (l3) are independent, the distributions 

in figures 1(a) and (b) apply for either a constant or variable wall tem- 
peratxare and for all Prandtl numbers. Increasing coolant flow (j^wl 

creasing) reduces the velocity at a given value of n; this amounts to a 
reduction in local shear stress. The distribution for f = -1.2 in 

V 

figure l(a) srlelds very small values for the velocity gradients (e.g., 
f"(0) = 0.0033575). When f^ » -1.23849, the velocity gradient at the 

wall becomes zero, according to reference 20. 

The S-shai>e of the velocity profiles, undesirable with respect to 
stability, is quite marked in figure l(a) for all nonzero values of f^. 

Inspection of figure l(b) for stagnation flow reveals that even for large 
valiaes of blowing there is no inflection point inside the boundary layer. 
This stabilizing influence of the pressure gradient is discussed in con- 
siderable detail in reference 22. 

Figiares l(c) and (d) show the effects of variable property values 
on the velocity distribution when the wall temperature is constant. The 
different values of T ^^T^ do not appreciably affect the shape when 

Eu = Oj for steignation flow (Eu = l), however, there is a marked effect 
of T ^T^ on the distribution, ivhich actually results in velocity over- 

^ shoot > 1^ for T j^T^ < 1. For the case of a favorable pressure gra- 

00 

dlent and a heated wall (Eu = 1, T^^yT^ l) ^ "the density is lowered so 
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that the flow iasidfi the boundary layer is accelerated more than the ex- 
ternal flcrw is accelerated, by the press\zre forces; overshoot results. 

For large tenrperature differences, the main stability criterion is 

^ (p • Figure 1(c) shows there are no inflection points in- 
side the boundary layer regardless of temperature ratio. 


Figures l(e) and (f) illustrate the combined effects of a porous wall 
and variable property values for a constant wall temperature. The veloc- 
ity gradients (and, hence, shear stress) are decreased by increasing 


Kw • 


Comparison of the curves for Eu = 0 with those for Eu = 1 shows 
that in all cases the pressure gradient tends to reduce the boundary- 
layer thickness and to increase the shear stress. 


Figure 2 contains plots of the temperature-difference ratio 
(T - Ty)/Too - T^^) = 0 against the dinffinsionless boundary -layer coordi- 
nate T|. Figure 2(a) shows these temperature profiles for various coolant- 
flow rates throu^ a porous flat plate (Eu =0) for constant property val- 
ues and zero ten^^erat-ure gradient at the wall, that is, for » 1 

and ®.jj.(0) 0* The values of n corresponding to the various coolant- 

flow rates are indicated in the figure. It can be seen in figure 2(a) 
that increasing the flow rate through the porous plate (i.e., increasing 
|f-;.r|) flattens the temperatvcre profile appreciably and increases the 

thickness of the temperature boundaiy layer. Similar profiles for stag- 
nation flow (Eu = 1) are shown in figure 2(b). Here, the effect of an 
increase in |f.^j.j is not nearly as pronoimced as for the flat plate (fig. 

2(a)). A conparison of figures 2(a) and (b) reveals the effects of Euler 
number on the temperature profiles for the case imder consideration. In- 
creasing the Euler number increases the temperature gradient in the bound- 
ary layer for each value of f^, resulting in decreased thickness of the 

temperature boundary layer. 

Figures 2(c) and (d) present temperature profiles for Eu = 0 and 
Eu = 1, respectively, for various coolant-flow rates and both constant 
property values and wall temperature, that is, for ^ n = 0. 

Curves representing profiles for Prandtl numbers of 0.7 (ref. 6) and 1 and 
2 (ref. 10) are also given for Eu = 0. Results for Pr = 1.0 with stag- 
nation flow (Eu = 1, fig. 2(d)) are taken from reference 2. 

For flat-plate flow (fig. 2(c)) the different effect of Prandtl num- 
ber for f = 0 and -1 is noteworthy. For the impermeable plate, in- 

IT • 

creasing Pr steepens the temperature distribution. For the strongly 
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cooled plate = -l), increasing Pr reduces the temperature near the 

wall. It should also he noted that in figure 2(d) the curves for Pr = 1 
with f^ = -3.1905 and -4.3346 resulted in a zero temperatxire gradient 

at the wall (d ’ (o) = O) . 

Curves for the flat-plate (Eu = O) and stagnation flows (Eu = l) are 
plotted in figures 2(e) and (f), respectively, for constant property val- 
ues, a Prandtl number of 0.7, and a linear wall-temperature distribution 
^ (n = l) . Ihe effects of increased coolant flow and of increased Euler 

H number are similar to those for a Prandtl number of 0.7 and a constant 

wall temperature, as shown in figures 2(c) and (d) . Compaidson of fig- 
ures 2(c) and (e), or figures 2(d) and (f), or "toth shovrs the influence 
of increasing the wall- temperature -gradient parameter n from 0 to 1. 

An increase in thfe temperature gradients in the boundary layer can be 
observed as n increases from 0 to 1. 

From the preceding discussion of figure 2, the following general 
trends can be noted: An increase in jf^j results in a decrease in the 

temperature gradients in the boundary layer for all values of n; an in- 
crease in n or in Eu, however, results in an increase in these tempera- 
ture gradients. These increases due to wall-temperature gradient are 
similar to those encountered in the velocity boundary layer due to main- 
stream velocity gradient (see fig. l) . The velocity boundary layer is 
affected by velocity gradients in the mainstream (outer edge of the bound- 
axy layer); the temperature boundary layer is Influenced not only by the 
velocity gradient but also by the temperature gradient along the wall 
(inner edge of the boundary layer) . 

Temperature profiles are given in figure 3 fa) for an impermeable flat 
plate (f^ = 0, Eu = 0) of constant temperature (n = O) for various stream- 

to-wall temperature ratios (large temperature differences) . Corresponding 
profiles for stagnation flow are shown in figure 3(b). Both figures show 
that Increases in the stream-to-wall hemperature ratio T reduce the 

temperature -gradient parameter 9 ’ in the boundary layer but have sU^t 
effect on the gradient at the wall. Euler number effects are similar to 
those noted in the previous figures. 

The effect of variations in fluid properties and of flow throu^ a 
porous wall are Illustrated for a flat plate and stagnation flow in fig- 
ures 3(c) and (d), respectively. For the flat plate, figure 3(c) shows 
that, over the range of t] plotted, the influence of Increased coolant 
flow through the plate exceeds that of increased temperature ratio; when 
Eu = 1, this situation is valid only for tj < about 2.2. 
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Examples of speciflc-wei^t-flow distri’butions for vaidalDle property 
values are given in figure 4 for constant wall temperatiore and stagnation 
flowj figure 4(a) is for an impermeable wall (f^ = O), and figure 4(b) is 
for a large blowing rate (f^ = -l.O) . For these cases, overshoot 

(—2^ > ij is encountered when > 1- The maximum overshoot occurs 

VPoo^oo ' 

for the permeable wall with maximum tenrperature ratio. 


fH 

a 

Friction and Heat Transfer Pj 

The dimensionless skin-friction parameter / 2 is sbo'vm as 

a function of pres sure -gradient paran^ter in figures 5 and 6 for imperme- 
able and permeable walls- It shorild be recalled that for small tempera- 
ture differences the velocity distributions are independent of n; there- 
fore, the same friction parameter will obtain for an wall-temperature 
variations. For large ten 5 >erature differences (T^^^T^ / however, the 

results apply only for constant wall tenperature. 

The influence of Eu on the friction is q\d.te marked, especially 
for Tj^T^ <1. For example, when T^T^ = l/4 and f ^ = 0 (fig. 5), 

an eightfold increase in the friction results when Eu is increased from 
0 to 1. Although T^T^ does not change the friction much for the flat 
impermeable plate, large changes are wrou^t for stagnation flow (Eu = l) . 
(similar results for Eu = 1 are ^ven in ref. 23. Son^ of the Ifech niun- 
ber aspects in ref. 23 are discussed in ref. 24.) The curves for differ- 
ent temperature ratios cross at Euler number near zero. The results in 
figure 5 are conpared with an approximate and simpler solution in refer- 
ence 25 (p. 47). The approximation is poorer for Eu < 0 than for 
Eu> 0- Larger adverse pressure gradients can be tolerated for strong 
T 

cooling (i.e., — > l) before the boiindary layer separates from the wall; 
separation, therefore, is delayed. 

Coolant-flow emission greatly reduces the friction for both small and 
large temperature differences (fig. 6) . In contrast to the flat imperme- 
able wall, a change in T^T^ from 1 to 4 reduces the friction by about 

one -third when f = -1.0. 

w 

Dimensionless local heat-transfer results are presented in figure 7 
for small temperature differences and variable wall temperature. A par- 
ticular wall -temperatiore variation results in a zero temperature gradient 
at the wall for each Euler number and coolant flow in figure 7. These 
values of n are given by the intercepts of the variooos cooives with the 
horizontal axis and moore accurately in table II . 
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For fixed values of Eu and f^, increases in n yield increases 

in ttie local convective heat-transfer parameter a,s a result of the in- 
creased gradients in the temperatirre profile. Figures 7(a) and (c) show 
that, for a linear wall-temperature gradient (i.e., n = l), a coolant 
flow represented approximately by f^ = -0.5 is required to obtain about 
the same heat-transfer coefficient as for a solid wall with constant 
temperature . 

As Eu increases from 0 to 1, Nu/V^ increases considerably -vdien 
n and f^ are fixed. Exceptions can be noted, however. The curves 

for stagnation flow (Eu = l) and a cooled wall (f^ = -0.5 and -l.O) are 
essentially the same as the corresponding curves for Eu = 0.5. This 
similarity emphasizes that the primary pressixre -gradient effects occur 
as Eu changes from 0 to 0.5. 

In flginre 8 the dimensionless local heat-transfer parameter is 
plotted against the Euler nttmber for large temperatTxre differences and 
a constant -temperature wall. The ciarves group themselves according to 
the coolant -flow rates considered, maximum coolant flow yielding minimum 
heat transfer. These ctrrves indicate the inadvisability of Increasing 
the coolant flow much more than that represented by a value of f„ of 
-1.0. Reductions in heat transfer accompany increases in temperature 
ratio T these reductions are especially marked when the wall is 
porous. Comparison of figures 5 and 8 shows that both friction and heat 
transfer are affected more by ten5>erature ratio for stagnation flow than 
for flat -plate flow. 

Heat-transfer results are shown in figure 9 for fluids with different 
Prandtl numbers and small temperature differences flowing over permeable 
and impermeable, isothermal flat plates (n = 0 = Eu) . Prandtl number 
ranges for several flxiids are Indicated along the abscissa. The results 
are converted from those given by Mckley and associates (ref. lO), by 
Eckert (refs. 2 and 8), and by Pohlhausen (ref. 7). There is good agree- 
ment between the present resiilts, those of reference 2, and the compar- 
able results of MLckley for f^ < 0. Common values are also obtained by 

the investigators for f ^ = 0 when calculations were made for identical 
Prandtl numbers . 

The formula of Pohlhausen 

= 0.332 (Pr)^/^ 

/\/Be 

is indicated on figure 9 and shows good agreement with the results when 
f^ = 0. When f^ <0, the results do not follow a simple power law. For 

f ^ — 0 the heat transfer is increased; for f ^ < 0 the heat transfer is 
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decreased, as already shown in figures 1 ' and 8 for Er = 0.7. Reference 
6 notes that the heat transfer is the same for f y = -0.5 when Pr is 

0. 7 and 1-0. From figure 9 it may he concluded that this similarity holds 
quite irell for 0.6 < Pr <1.5. Qhe marked reduction in heat transfer 
with coolant emission is especia3dy striking at the higher Prandtl num- 
bers. At a Prandtl number of 5.0, Wu/VRe is reduced by the ratio 
1/8250 when f^ changes from 0 to -1. This behavior may be due to the 
increased heat capacity of fluids with larger Prandtl numbers. In con- 
trast to results for the impermeable plate, the heat transfer decreases 
with increasing Prandtl number when f^ < 0. (The case f^ = -0.5 is 

exceptional over part of the range.) 

SUMMARY OF RESULTS 

Solutions of the land nar-boundary-layer equations were summarized 
for ranges of mainstream pressure gradient and rates of coolant emission 
through a porous wall- For small differences between the wall and stream 
temperatvure, the wall temperature was allowed to be either constant or 
variable j for large differences between stream and wall ten^erature only 
a constant wall ten^jerature was considered. Solutions obtained herein 
were restricted to a Prandtl number at the rail of 0.7 and negligible 
temperature changes caused by con^ression and frictional heating compared 
with the difference between the wall and mainstream temperatures. Dimen- 
sionless forms of heat-transfer and friction parameters and boundary-layer 
thicknesses are tab\alated. 

The results of this analytical study are summarized as follows: 

1. The velocity and temperature distributions near the 'vrall indicated 
reduced gradients for increasing coolant emission and for increasing ratio 
of stream-to-wall ten^rerature . These gradients increased for increasing 
mainstream pressure gradients. For small temperature differences between 
wall and stream, the velocity greidients are unaffected but the teraperatvcre 
gradients are increased by increasing wall-temperature gradient. 

2. For stagnation flow over an impermeable wall, the boundary-layer 
velocity exceeded the mainstream velocity for stream-to-wall temperature 
ratios less than 1. For stream-to-wall temperature ratios greater than 

1, the specific wel^t flow in the boundary layer exceeded that in the 
mainstream for stagnation flow over either a permeable or an Inpermeable 
wall. 


3. Stream-to-wall temperat-ure ratios above 1 and coolant emission 
diminished the friction. This effect was more marked for stagnation flow 
than for flat-plate flow. 

4. Increased convective heat transfer accompanied larger wall tem- 
perature gradients. The heat transfer was reduced, however, by coolant 
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emission through the wall. These results are applicable for flat-plate 
and stagnation flows with small temperature differences (constant property 
values). For large temperature differences and a constant wall ten^jera- 
ture, the heat-transfer parameter was reduced for increasing stream-to- 
wall ten^jerature ratios. This effect is especially pronounced -vdoen the 
wall is porous. 

5. Results were found in the literature for flow of different flxoids 
with small temperature differences along permeable and impermeable iso- 
thermal flat plates. These resTilts indicated much greater reductions in 
heat transfer due to coolant flow for Erandtl number of the order of 5 
than for Erandtl number of the order of 0.7. In contrast to results for 
the impermeable plate, the heat transfer decreased with increasing Erandtl 
number and transpiration cooling. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 11, 1955 
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TABLE !• - lAMim-BCURDARY-IAm SOIOTIOHS 


Investigator 

Ref. 

Type of 
solvcbioD 

Pressure 

gradient 

Prandtl 

number 

Wall 

txire 

Type 

wall 

Property 

variation 

Pohlhausen 

7 

Exact 

Ro 

0.7<Pr< 15 

Constant 

lnpermaable 

Cons'tant 

Pallmer and Skan 

15 

Exact 

Yes 


— 

Impermeable 

Constant 

Page and PeiUmar 

16 

Exact 

Yes 

0.77 

Verialle 

Impermeable 

Cons'tant 

Eckert and Drewlte 

1 

Exact 

Yes 

0.7 

Constant 

Impermeable 

Cons'tant 

Eckart 

8 

Exact and approx. 

Yes 

0.6< Pr< 10 

Constant 

lapermeable 

Constant 

achllchtlng and BioBsmann 

21 

Exact 

Yes 

— 

— 

Permeable 

Constant 

Bchiih 

17 

Exact and approx. 

Yes 

0.7 

Variaile 

Impenneabla 

Cons'tant 

Eckert 

2 

Exact 

Yes 

1.0 

Constant 

Permeable 

Cons'tant 

Brcfwn and Donoughe 

14 

Exact 

Yes 

0.6<Pr<l.l 

Constant 

Impermeable 

Constant 

Brown 

3 

Exact 

Yes 

0.7 

Constant 

Permeable 

Variable 

Brcvn and Dononghe 

4 

Exact 

Yes 

0.7 

Cons'tant 

Permeable 

Variable 

Btanlf orth 

12 

Approx. 

Yes 

O'. 71 

Constant 

Permeable 

Cons'tant 

Brown and Liv ingood 

5 

Exact 

Yes 

0.7 

Constant 

Permasble. 

Variable 

Levy 

15^ 

Exact 

Yes 

U.YdiriiJU 

variance 

impermeable 

uonsnanx 

Eckert and Idvingood 

n 

Approx. 

Yes 

0.7 

Cons'tant 

Permeable 

Variable 

Eckert and Idvingood 

13 

Approx. 

Yes 

0.7 

Cons'tant 

Permeable 

Variable 

Tifford and Chti? 

9 

Exact 

Yes 

0.5<Pr< 10 

Variable 

Impermeable 


Emnona and lei^ 

20 

Exact 

Ho 

— 


Permeable 

Constant 

Thurston®' 

25 

Approx. 

Yea 

0.7 &1.0 

Constant 

Impermeable 

Variable 

Donongiie and Idvingood 

6 

Exact 

Yas 

0.7 

Variable 

Parmoable 

Cons'tant 

Mlckley 

10 

Exact 

Ho 

0.6<Pt55.0 

Cons'tant 

Perrmeabla 

Cons'tant 

Levy 

23 

Exact 

Yes 

0.7 &1.0 

Cons'tant 

Impermeable 

Variable 


^Siumnnry rejKirt. 
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TABI£ II. - SOMMARI OF HEAT-TRAHHFER AMD FKECTIOH PARAI.EEERS AHD 
BOUHDARX-LAYER THICKIlESfflS 


(a) Small temperature differences s? » 1 







X 

X 

0.626 

2.737 

.692 

2.505 

.720 

2.415 

.752 

2.307 

.773 

2.233 

.801 

2.128 

1.861 

3.417 

.835 

1.959 

.578 

1.562 

.458 

1.363 

1.6075 

2.403 

.792 

1.407 

.620 

1.186 

.518 

1.0515 

1.409 

2.0175 

.9185 

1.440 

.708 

1.187 

.586 

1.036 

.506 

.936 

0.807 

3.677 

1.921 

4.134 

.974 

2.6495 

.623 

2 ..047 

.472 

1.763 

1.743 

2.706 

1.651 

2.593 

.994 

1.776 

.738 

1.445 

.599 

1.258 

1.5535 

2.288 

.918 

1.530 

.7265 

1.292 

.609 

1.143 

1.072 

5.997 

1.909 

5.860 

1.150 

4.414 

,650 

3.338 

.465 

2.871 

1.241 

3.317 

1.280 

2.806 

1.905 

3.078 

1.270 

2.274 

.889 

1.783 

.700 

1.524 

1.731 

1.208 

.910 

2.622 

1.995 

1.627 

.740 

1.411 
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TABIE U. - Concluded. SDl.lMAEy OF EEAT-TRAHSFER AHD FRICTION PARAMETERS AND 

BOCHDARY-LAm THICKNESSES 

(Ij) Large teinperature differences with constant wall temperatva-e (n » 0) 



"C. 

Tv 

Eu 

Hu 

VRe 

6' 

w 

Cf.wVS 

2 

f" 

w 

8*4^ 

X 

X 

8g -VEe 
X 

8^ VSe 

X 

Additional 

infonnatlon 

Refer- 

ence 

Page 

0 

0.25 

0 

0.2884 

0.3556 

2.031 

0.182 

0.246 

1.482 

1 

1 

62 



.50 

.4801 

1.9299 

1.033 

-.031 

.270 

.919 



64 



1.00 

.5812 

2.7842 

.820 

-.064 

.246 

.764 



65 


0.5 

-0.06 

0.2062 

0 

3.043 

0.442 

0.348 

2.215 


L 

57 



-.04 

.2554 

.1735 

2.309 

.406 

.420 

1.846 



58 



0 

.2900 

.3462 

1.898 

.347 

.457 

1.645 



59 



.50 

.4412 

1.2754 

.980 

.116 

.460 

1.109 



60 



1.00 

.5298 

1.8000 

.768 

.065 

.415 

.929 



61 


2 

-0.1178 

0.1890 

0 

4.582 

1.663 

1.076 

3.721 

k 

L 

16 



-.09 

.2522 

.1634 

2.430 

1.501 

1.408 

2.928 



17 



-.05 

.2756 

.2434 

1.882 

1.378 

1.478 

2.713 



18 



0 

.2944 

.3125 

1.537 

1.271 

1.495 

2.547 



19 



.50 

.4020 

.6836 

.622 

.838 

1.372 

1.904 

5 

18 



1.00 

.4760 

.9090 

.433 

.693 

1.215 

1.611 

5 

19 


4 

-0.1351 

0.1794 

0 

6.950 

3.100 

1.834 

5.596 



22 



-.09 

.2642 

.1954 

2.297 

2.719 

2.595 

4.053 



23 



-.05 

.2790 

.2397 

1.810 

2.582 

2.651 

3.863 



24 



0 

.2952 

.2874 

1.438 

2.455 

2.663 

3.678 



25 



.50 

.3940 

.5541 

.272 

1.676 

2.399 

2.790 

5 

24 



1.00 

.4662 

.7220 

.081 

1.403 

2.113 

2.370 

1 

> 

25 

-o.s 

2 

0 

0.1602 

0.1476 

2.381 

1.605 

1.778 

3.456 

4 

32 



.50 

.2304 

.4770 

.772 

1.044 

1.771 

2.476 

5 

20 



1.00 

.2560 

.6415 

.507 

.880 

1.638 

2.151 

5 

21 


4 

-0.0644 

0.0796 

0 

7.219 

3.485 

2.620 

7.329 

4 

35 



0 

.1506 

.1263 

2.460 

3.100 

3.237 

5.012 

4 


36 



.50 

.2088 

.3448 

.318 

2.155 

3.142 

3.670 

C 

s 

> 

26 



1.00 

.2252 

.4465 

0 

1.867 

2.932 

3.253 

C 


27 

-1.0 

2 

0 

0.0406 

0.0242 

4.931 

2.114 

2.167 

5.941 

4 


45 



.05 

.0692 

.0892 

2.985 

1.911 

2.299 

4.625 



46 



.15 

.0686 

.1678 

1.989 

1.704 

2.343 

3.959 



47 



.50 

.1062 

.3228 

.965 

1.310 

2.298 

3.236 

C 

) 

22 



1.00 

.1052 

.4445 

.608 

1.148 

2.228 

2.891 

5 

23 


4 

0 

0.0262 

0.0125 

6.409 

4.160 

4.002 

9.007 

4 


51 



.05 

.0510 

.0542 

3.405 

3.881 

4.247 

6.995 



52 . 



.15 

.0656 

.1030 

2.040 

3.620 

4.309 

6.054 



53 



. .50 

.0780 

.1959 

.402 

2.882 

4.256 

4.954 

C 


28 



1.00 

.0726 

.2611 

-.117 

2.602 

4.147 

4.448 

c 
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(c) Plat-plate flow for constant wall (d) Stagnation flow for ccnatant waill 

temperature. Eu, 0; n, 0. temperature. Eu, 1.0; n> 0. 



0 1234560125456 

Nondlraensional boundary-layer coordinate, 77 


(e) Flat-plate flow for linear wall-temperature (f) Stagnation flow for linear wall-temperature 
variation. Eu, 0; n, 1.0; Pr, 0.7. variation. Eu, 1.0; n, 1.0; Pr, 0.7. 

Figure 2. - Temperature distributions in laminar boundary layer with small 
temperature differences » l). 
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Figure 6. - Skin- friction parameter for laminar lioundary layer for 
permeable and impermeable Trails at constant temperature. Wall- 
temperature-gradient parameter, 0. 
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Wall- t*op«i«tur*-gradaent paro«ttr, n 
(a) Eulor tnuibar, 0. (b) Hular nuobsr, O.E. 

FIsum 7. - Baat-trauafer para»etar for ooa»tant-pi>oparty laalnar boundary layer for 
Pmndtl ntaaior, 0,7. 
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Euler xmmber, Eu 


Pigure 8. - Heat-transfer paraseter for laxalnar boundary layer for 
permeable and Impermeable -valla at constant temperature. ¥all- 
temperature-sradient parameter, 0; Frandtl nuab^, 0.7. 
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Figure 9. - Heat-transfer parameter for constant 
property laminar boundary layer for different 
fluids. Constant wall temperature, wall- 
temparature-gradlent parameter, 0; Euler num- 
ber, 0. 
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